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We calculate the heat flux and electron-phonon thermal conductance in a disordered graphene 
sheet, going beyond a Fermi's Golden rule approach to fully account for the modification of the 
electron-phonon interaction by disorder. Using the Keldysh technique combined with standard 
impurity averaging methods in the regime kpl 3> 1 (where kp is the Fermi wavevector, I the 
mean free path), we consider both scalar potential (i.e. deformation potential) and vector potential 
couplings between electrons and phonons. We also consider the effects of electronic screening at 
the Thomas-Fermi level. We find that the temperature dependence of the heat fiux and thermal 
conductance is sensitive to the presence of disorder and screening, and reflects the underlying chiral 
nature of electrons in graphene and the corresponding modification of their diffusive behaviour. In 
the case of weak screening, disorder enhances the low-temperature heat flux over the clean system 
(changing the associated power law from to T'^), and the deformation potential dominates. For 
strong screening, both the deformation potential and vector potential couplings make comparable 
contributions, and the low-temperature heat flux obeys a power law. 

PACS numbers: 65.80.Ck, 72.10.Di, 44.10.-|-i 



I. INTRODUCTION 

The potential to exploit the exceptional thermal prop- 
erties of graphene in applications has recently gener- 
ated considerable activitjfflEl possible applications in- 
clude sensitive bolometry and calorimetry for detecting 
infrared and THz radiation. Such detectors would ul- 
timately be based on the simple heating of electrons in 
a graphene sheet by the absorption of incident photons. 
Ideally, the detector electrons would be thermally decou- 
pled from their surroundings, thus allowing any heating 
produced by the incident radiation to be long-lived. In 
this respect, graphene provides a potential advantage: its 
low electron density and relatively weak electron-phonon 
coupling implies that the expected low-temperat ure ther- 
mal decoupling between electrons and the latticd^ could 
occur over a muc h w ider temperature range than in a 
conventional metal^. Further, one can effectively sup- 
press the thermal link between graphene electrons and 
electrons in the contacts by employing superconducting 
leadiP. 

Given the above, it is crucial to develop a rigorous 
and quantitative understanding of the thermal link be- 
tween electrons and phonons in graphene at low temper- 
atures. Several recent theoretical works have addressed 
this problem in the ca se of clean graphene (i.e. no elec- 
tronic disorder jP^^SESl, At low temperatures, the heat 
flux (per volume) between electrons and longitudinal 
acoustic phonons takes the general form: 

P(r„ Tph) = F(Te) - F(rph) = S(T<f - Tp^J, (1) 

where F{T) is called the energy control function, S is 
a coupling constant, and To and Tph are the tempera- 
tures of the electrons and latt ice (i.e. phonons) respec- 
tively. Previous workJ^l^Mini that 5 — A in the 
low-temperature limit, assuming an unscreened deforma- 
tion potential electron-phonon coupling. This is identical 



to what would be expected for a clean conventional 2D 
metaM 

In this work, we now ask how the above result is 
modified in the presence of electronic disorder. While 
great experimental progress has been made in reducing 
disorder effects in graphen^li', the devices studied for 
bolometric applications in Refs. 12)31 are sitting on a sil- 
icon substrate and have mean free paths that are 100 
nm or less. In conventional metals, electronic disorder 
can strongly affect the electron-phonon coupling at tem- 
peratures low enough that the wavelength of a thermal 
phonon is comparable to (or longer than) the electronic 
mean free path^^ This defines a characteristic tem- 
perature scale Tdis below which disorder effects are im- 
portant, 

fceTdis = hs/l, (2) 

where s is the speed of sound and I is the electronic mean 
free path. As discussed in Refs. [T^HTSl (and below), the 
effects of disorder are subtle: depending on the nature 
of the disorder and the electronic system, the power law 
6 in Eq. ^ can either be enhanced by disorder or be 
suppressed^21tl5| 

Here, we study how the additional complexity arising 
from the unique electronic properties of graphene modify 
the interplay of disorder and the electron-phonon inter- 
action. The principle new ingredients arise from the ef- 
fective chiral nature of carriers in graphene, which both 
modifies electronic diffusion, and allows for a new kind of 
effective vector-potential electron-phonon coupling^^. 
For simplicity, we will focus on disorder originating from 
charges in the substrate below the graphene fiake, and 
thus take the disorder potential to preserve the symme- 
tries (valley and sublattice) of the low-energy graphene 
HamiltoniarfiSlini, The impurity potential is thus also 
taken to be static, i.e. it does not move with the graphene 
sheet. We also consider the case where the graphene 
fiake has been doped sufficiently that kpl ^ 1 {kp is the 
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Deformation potential 


Vector potential 


T < Tdis 
weak screening strong 


Tdis <T <Tbg 
' screening weak screening strong screening 


T < Tdie Tdis <T <Tbg 
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TABLE I: Energy control function of deformation potential and vector potential in graphene below the Bloch Griineisen 
temperature Tbg = 2hskF /k^- gi and g2 are the deformation potential and vector potential coupling constant respectively 
(c.f. Eqs. ([6|). Ef is the Fermi energy with respect to the Dirac point, vf the Fermi velocity, I the electronic mean free path, 
and pM is the mass density of graphene per area, gxp is the Thomas-Fermi wavevector (as defined in Ref. |23p. and ({n) is the 
zeta function. 



Fermi wavevector) , meaning that the standar d imp urity- 
averaged perturbation theory is appropriate! ^^ * ^^ ! Com- 
bining this approach vifith the Keldysh technique then 
allows us to rigorously address the electron-phonon in- 
teraction in the presence of disorder, in a manner anal- 
ogous to the classic works looking at this physics in a 
conventional disordered metaP^HEl^ stress that prop- 
erly addressing disorder effects involves going beyond the 
sort of Golden Rule calculation used to address the clean 
cas^. 

Note that since we work in the regime kpl ^ 1, 
the temperature scale Tdis below which disorder effects 
emerge will necessarily be well below the Bloch Griineisen 



temperature feBlBG = 2hskF', we will thus explicitly fo- 
cus on temperatures T < Tbg- For a typical doped 
graphene electron density 10^^/cm^, Tbg ~ 70K. In 
contrast, for a typical mean free path of lOOnm and 
a graphene acoustic phonon velocity of ~ 2 x lO^m/s, 
Tdis ^ lOK. We note that a recent study examined disor- 
der effects on electron-phonon interactions above the BG 
temperature^. Somewhat surprisingly, the expression de- 
rived in that work for impurity-assisted electron-phonon 
cooling for T ^ Tbg is exactly half of our expression for 
e-phonon cooling based on the deformation potential at 
T < Tdis < Tbg (cf. Table |l] and Eq. ([38])). 



Our main results for the energy control function F(T) 
(c.f. Eq. ([T])) for graphene at low-temperatures are sum- 
marized in Table IH We consider the contribution to the 
electron-phonon heat flux arising from both the standard 
deformation potential coupling (DP), as well as from the 
effective vector potential coupling (VP). In the absence 
of disorder (or at temperatures well above Tdis but be- 
low Tbg), &nd in the absence of electronic screening, one 
finds that both these mechanisms contribute indepen- 
dently and in a similar manner: the respective heat fluxes 
are each described by Eq. ([T]) with 6 = 4: (in agreement 
with Refs.mni). 

In contrast, for T < Tdis, the two coupling mecha- 
nisms are affected oppositely by the electronic disorder. 
We find that the heat flux associated with the VP cou- 
pling is suppressed by disorder: still neglecting screening, 
it is described now by Eq. (!]) with an enhanced power- 
law of (5 = 5. Heuristically, this is attributed to the 
disorder-broadening of the graphene energy levels. The 
effect of disorder on the heat flux associated with the DP 
coupling for T < Tdis is the opposite from the above: 
it is enhanced. It is described by a reduced power-law 
6 = 3 (again, no screening), and will thus dominate the 
VP at low temperatures. On a heuristic level, this en- 
hancement is due to the diffusive charge dynamics, which 



effectively increases the time an electron interacts with a 
given phonon (i.e. this becomes the time to diffuse across 
a phonon wavelength, as oppose the time needed to bal- 
listically traverse this distance). The absence of any dif- 
fusive enhancement of the VP coupling is (as we will 
show) a direct consequence of the non-conservation of 
pseudospin. 

We also consider how including screening changes the 
above results; the importance of screening the e-ph inter- 
action has been the subject of several recent studieJ^SHUl, 
As discussed extensively in Refs. I23|25|, the DP coupling 
is expected to be screened, whereas the VP coupling is 
expected to be unscreened, as it induces no net electronic 
charge (i.e. the effective vector potential generated by a 
phonon field has opposite sign in the two graphene val- 
leys). As a result (see Table I), even without disorder 
(i.e. Tdis < T < Tbg), the DP and VP heat fluxes are 
not equivalent in the limit of strong screening: the VP 
power-law remains 6 = A, where the DP power law is 
increased to 6 = 6. Similarly, at low temperatures where 
disorder effects matter, the VP power-law is unchanged, 
but the DP power law becomes T^. 

Note that our results suggest that even though the 
bare VP coupling strength 52 is believed to be about 
an order-of-magnitude smaller than the bare DP cou- 
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pling g^^, if screening is strong, its contribution to the 
heat flux could be comparable to or even large than that 
from the DP coupling. This is despite the relative en- 
hancement of the DP coupling over the VP coupling by 
disorder. Further, our results suggest that the electron- 
phonon heat flux could be a means for empirically de- 
termining if screening is important. In particular, the 
only way to obtain a power-law is via an unscreened 
DP coupling in the diffusive limit. We note that mea- 
surements of the phonon contribution to the resitivity in 
a clean graphene sheet (as recently measurecP^ cannot 
directly resolve this issue, as both DP and VP couplings 
contribute a dependenc^^U^ 

This paper is organized as follows. In Sec. |llj we 
present our model and an outline of the calculational 
method. This includes a brief derivation of the kinetic 
equation of electrons in graphene in the Keldysh formal- 



ism (Sec. II B ), as well as a derivation and discussion 
of the diffusion propagator in graphene and the resulting 
diffusive renormalization of the electron-phonon vertex 
(Sec. II C). We present the main results in Sec. |III[ i.e. 
the heat flux due to electron phonon interaction in both 
weak and strong screening case. Finally, we briefly sum- 



marize the paper in Sec. IV 



II. MODEL AND CALCULATION 

A. Model Hamiltonians of impurity and 
electron-phonon scattering 

1. Electrons m disordered graphene 

The low-energy electronic degrees of freedom are de- 
scribed by a massless Dirac Hamilton ian. Focus ing on a 
single valley (the K+ valley), one has' " ' ^^ ' ^^ !^ 

%= f (fr^\r){-ivF&jdj + U{r)i)<i{Y) (3) 



is a spinor field operator describing the am- 



where vp — lO^m/s is the Fermi velocity, \I'(r) 
( ^A(r) ' 
V ^B(r) 

plitudc of electrons on the two sublattices, Uj {j — x, y) 
are Pauli matrices, and U{r) is the disorder potential; 
we also set h = 1 throughout unless otherwise indicated. 
We do not include an index for spin or valley, as for the 
physics we consider, each spin and valley contributes in 
an equal and independent fashion. 

As mentioned, we focus on a smooth disorder potential 
originating with impurities in the substrate. We thus 
treat the impurity potential U (r) as a scalar potential 
with respect to both the valley degree of free dom a nd 
the sublattice degree of freedom (i.e. pseudospin^P^Ml, jn 
the standard way, U will be treated as delta-correlated 
Gaussian disorder, with zero-mean and correlator: 



The corresponding scattering rate is 1/t = vp/l = t^vw, 
where v = kp/i'KVF is the density of states at the Fermi 
energy per spin per valley. 



2. Electron-phonon interaction 

The electron-phonon interaction in graphene has been 
studied extensively in Ref. I23I25I For suspended 
graphene, there are both in-plane phonon modes and fiex- 
ural phonon modes (out-of-plane). In this work, we con- 
sider graphene on a substrate (as in recent experiments 
probing thermal propertieJ^El), such that flexural motion 
is suppressed; we thus only focus on in-plane motion. 
Further, at low to moderate temperatures, the optical 
modes are barely excited and the dominant modes par- 
ticipating in cooling of hot electrons are acoustic modes; 
we thus focus exclusively on the coupling to these modes. 

Due to the Dirac Hamiltonian of electrons in graphene, 
there are two distinct electron-phonon coupling mecha- 



The first is a standard deformation potential 
coupling, which corresponds to a local dilation of the lat- 
tice. In the Dirac theory, it appears as a scalar potential 
(with respect to pseudospin). The second mechanism is 
an effective gauge-field coupling or vector potential cou- 
pling. This corresponds to the change in hopping ma- 
trix elements accompanying a pure shear deformation, 
and enters the Dirac theory the same way as an exter- 
nal gauge field (the only proviso being that this phonon- 
induced vector potential is valley-odd, and hence does 
not break time- reversal symmetry). 

Letting 'S'k denote a momentum-space electronic field 
operator and 5^ q a phonon annihilation operator, the 
total interaction between electrons and acoustic phonons 
can be written in the general forrrP^ 

= E E *k+q A^'lq) *k (&,,q + hl_^ . (5) 
r}—l,t k.q 

Here T] = l,t denote longitudinal (LA) and transverse 
(TA) acoustic modes respectively. The 2x2 coupling 
matrices M'^{q) take the form: 



M*(q) = iqC, 



W2 



9i 
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where 



526 



(6a) 
(6b) 

(7) 



(C/(r)t/(r')) ^w6{r~r'). 



(4) 



Here pm is the mass density of the graphene sheet, is 
the phonon frequency for ij mode, is the angle of the 
phonon wavevector q with respect to the x axis (which is 
taken to be along the armchair direction of the graphene 
lattice) . For simplicity, we take the speed of sound to be 
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the same for LA and TA phonon and drop the superscript 
rj in the frequency from now on. gi (172) is the deforma- 
tion potential (vector potential) coupling constant. Pre- 
vious works have estimated gi ^ 20 — 30 eV and 32 ^ 1-5 
e\^, though we note that even the value of the deforma- 
tion potential coupling is subject to some debateP22Ml. 
Our theory is thus not tied to specific values of these pa- 
rameters, and we keep both the DP and VP couplings in 
our discussion. Note that transverse phonons induce only 
a pure shear deformation, and hence couple only through 
the vector potential. 



phonon 




(b) 

e,p, e,p/ e,p/ 

AAA*^ = ^^^^^^ ^ ^^^^^^^i 
e-|-a;,p + q e-|-a;,p-|-q e-|-a;,p-|-q 



B. Keldysh formalism of the kinetic equation of 
electrons in graphene 

Having established the basic electronic and electron- 
phonon Hamiltonians (c.f. Eq. ([s]) and ([5|), we now turn 
to our main goal of calculating the heat flux between elec- 
trons and phonons. We consider the standard situation 
where each subsytem is independently in thermal equi- 
librium at its own temperature (electrons at To, phonons 
at Tph). In the disorder-free case, this heat flux can be 
conveniently calculated by using Fermi's Golden rule to 
calculate electron-phonon scattering rate^^^^. Includ- 
ing disorder, we need a more general formalism, one that 
is capable of capturing the interference between electron 
phonon and electron impurity scatterings (i.e. the ver- 
tex correction of the electron-phonon vertices M''(q) by 
disorder) . To that end, we make use of the Keldysh 
techniqu d'^^ l ^^ l, coupled with standard disorder-averaged 
perturbation theory. Such an approach was used by 
Kechedzhi et al. to study conductance fluctuations (in 
the absence of any electron-phonon coupling)^. 

We start by noting that the heat flux of interest (i.e. 
energy lost/gained by the electrons) can be directly re- 
lated to the collision integral /o(£, p) appearing in a stan- 
dard Boltzman equation describing the dynamics of the 
electronic phase-space distribution function n[e, p;i]: 



dn[e, p] 



dt 



(8) 



-ph scatt 



The collision integral /o(£,p) tells us the rate of change 
of p; t] due to the emission and absorption of acoustic 
phonons. 

As electronic momentum relaxation is much faster than 
energy relaxation, to describe the latter process we can 
focus on times longer than the electron momentum re- 
laxation time; in addition, the distribution function will 
be sharply peaked on-shell (magnitude of momentum set 
by energy). The relevant kinetics can thus be described 
by an electron distribution function n{e;t) that depends 
only on energy, not on momentum. The kinetic equation 
for the electronic distribution function n(e;i) takes the 
form 



dn{e\ t) 
dt 



1 



dp 



/o(e,p)A(e,p), (9) 



FIG. 1: Color online. (a)Electron self-energy diagram. The 
wavy line represents phonon propagator, the solid line rep- 
resents electron propagator and the square block represents 
dressed electron phonon vertex by difluson. (b) Vertex correc- 
tion of the electron-phonon interaction by impurity scattering. 
Dashed line represents an impurity average. The square block 
represents dressed vertex and the dot represents bare vertex. 



where 



A{e, p) — — Im 



2(£+2V) 



(10) 



is the electron spectral function. 

The heat flux between electrons and lattice (for one 
valley and spin projection) is given by 



(11) 



The collision integral /o(£, p) is obtained in the stan- 
dard manner by calculating the electronic Keldysh self- 
energies S arising from the electron-phonon interaction, 
to first order. One finds the general relatiorP^ (see Ap- 
pendix IA|) : 



--Tr 



(12) 



The self-energies & {j = K, R and A) are 2x2 matrices 
(in pscudospin space), and are functions of both energy 
£ and momentum p; they include the effects of disor- 
der averaging. As we are considering a quasi-equilibrium 
situation where both phonons and electrons are individ- 



ually in thermal equilibrium, n(£, p; t) in Eq. ( 12 ) can 



be replaced by a Fermi distribution function at tempera- 
ture Tc, and the self-energies can be calculated assuming 
phonons are in thermal equilibrium at temperature Tph. 

In the regime of interest {kpl ^ 1), the dominant 
self-energy diagram describing the leading-order electron- 
phonon contribution to the kinetic equation (in the 
presence of disorder) is shown in Fig. [T] In this di- 
agram, the wavy line represents a phonon propagator 
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{oj , q; r]) , where /? = R,A,K denotes retarded, ad- 
vanced and Keldysh propagators. The retarded and ad- 
vanced phonon propagators appearing here are 



n«M(a;,q) 



2uj„ 



(13) 



with 



sq, while the Keldysh propagator is 



n^(L^, q) = (1 + 2N{uj, Tph)) (n«(c^, q) - W'iu;, q)) , 

(14) 

where N{Lj,Tp\^) is the Bose-Einstein distribution evalu- 
ated at T = Tph. 

The solid line in Fig.[T]represents an impurity-averaged 
electronic Green function G^{e,p); note that these are 
2x2 matrices in pseudospin space. The retarded and 
advanced components are given by: 



e± 



vpo ■ p 



(e± ^)2 - W^p2 
while the Keldysh electron Green function is 

G^(£, q) = (1 - 2n(e, T,)) (G«(e, q) - G^{e, q) 



(15) 



(16) 

where n(e, T) is now the Fermi-Dirac distribution func- 
tion evaluated at T = To and chemical potential Ep^ 
where Ep is the Fermi energy (measured from the Dirac 
point). 

Note that the electron Green function in Eq. (151 has 
an extremely simple form: it is just a free propaga- 
tor with the substitution e — ?> e ± i/2T, corresponding 
to disorder-induced broadening of energy levels. This 
broadening represents the first mechanism by which the 
electron-phonon heat flux will be modified due to dis- 
order; this broadening generally causes a suppression of 
the heat flux. The second key effect of disorder is via the 
vertex correction of the electron-phonon vertex appear- 
ing in the self-energy in Fig. [T] Each vertex describes 
the emission or absorption of a phonon; heuristically, 
the disorder-induced vertex correction corresponds to the 
modification of the amplitude of such a process due to 
the diffusive motion of electrons. The full details of the 
renormalization of the e-phonon vertices by disorder in 
the Keldysh formalism in normal metals are presented in 
Ref. |13j. The Keldysh structure and renormalization of 
the e-phonon vertices in graphene can be treated in a sim- 
ilar fashion. The only key difference comes from the 2x2 
matrix structure associated with pseudospin; as we will 
see, this leads to interesting new physical consequences. 

The electron phonon vertices M (the rj superscript is 
dropped here) in the Keldysh technique are represented 
by the form , where the upper index is for phonons 
and the lower electrons. The index 7, /i, v each has two 
components, cl and g, due to the two-component struc- 
ture of the electron and phonon fields in the Keldysh 
formalisrrP^. The vertices with different indices are 
rcnormalizcd differently by disorder in the Keldysh tech- 
nique as shown in Ref. 1131 Here, we only present the sim- 
plest case, the renormalization of the vertex M^^ Upon 



summation of all self-energy diagrams, one finds that the 
form of this renormalized vertex appears directly in the 
final expression for the collision integral, Eq. (21 1. The 



renormalization of vertex M^^ by disorder is depicted 
in Fig. [TJd. We focus on this specific vertex in the re- 
mainder of this subsection and drop the upper and lower 
index from now on. 

In general, we may describe the renormalization of a 
electron-phonon interaction vertex by 



Mdiff(q,w) = £'(q,w) o Afo(q) 



(17) 



where Mq (Mdiff) is the bare (renormalized) vertex, and 
-D(q, w) is a linear operator acting in the space of 2 x 2 
matrices. It represents the diffusion propagator for elec- 
trons in graphene, with the non-trivial matrix structure 
reflecting the fact that charge and pseudospin diffusion 
are linked together. It is convenient to write this ex- 
pression using a basis of Pauli matrices. Defining the 
four-component vectors m, mdiff via 



Mo = "iO ■ (i,'5'a;,'5'y,'5'z) , 

— "idiff • (1, (To;, (7y, (7z) 



Eq. (17) takes the form: 



mdiff (q, t^) = 2?(q, Lo) ■ ™o(q) 



(18a) 
(18b) 

(19) 



where I? is a 4 x 4 matrix. Using this representation, the 
summation of ladder diagrams depicted in Fig. [T] results 
in the form: 

2?(q,L^)-(l-7'(q,c^))^\ 



(2^ 



G^(k + q,cj) G^(k,0) 



ill 



'Ik- 



(20) 



Here, the indices a, /3 run from to 3, and repeated in- 
dices are to be summed over; we also use ct*^ denote the 
2x2 unit matrix and tr' (i = 1, 2, 3) are the Pauli matri- 
ces. V describes a single "rung" in a standard diffuson 
ladder. Its 4x4 matrix structure is now directly re- 
lated to the fact that each propagator carries an initial 
and final charge or pseudospin index. A similar structure 
is encountered when considering diffusive dynamics in a 
system with strong spin-orbit coupling or in 2D helical 



metals, as studied by Burkov et alP^ISS. 

We will discuss the properties and physics encoded 
in the matrix diffusion propagator T) in more detail in 
the next subsection. For now, we only show how it en- 
ters in the final expression for the collision integral (and 
hence the heat flux). One finds that due to the causality 
structure of Keldysh Green functions, the two vertices in 
Fig. [1^ cannot both be simultaneously dressed by impu- 
rity scattering. Summing up all the self energy diagrams 
in the Keldysh formalism, one finally obtains: 
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v=i,t L 



(27r)3 



M^'(q) n«(-L.;r,)G^(e + c^,p + q)M^ig(-c.,-q) + h.c, 



(21) 



where R{e,uj) is the expected combination of Bosc- 
Einstein and Fermi-Dirac functions appropriate for 
phonon emission and absorption processes: 

i?(£,w) = 7V(cj,Tph)n(e,Te)(l-n(£ + cj,re)) 

-(1 + N{io, Tph))(l - n{e, T,))n{e + T,) 
= {n{e, Te) - n{e + u, T,))[N{lo, Tph) - N{lo, T,)] . 

(22) 

The first term of R{e, oj) describes absorption of a phonon 
uj from energy state e to e + w and the second term de- 
scribes emission of a phonon lu from energy state e + uj to 
state e. As expected, R{e,uj) vanishes if Te = Tph. We 
note that apart from the matrix structure of the electron- 
phonon vertices and electron propagators, the expression 
for the coUision integral has the same form as that found 
for conventional diffusive metals^^ 13 15' Nonetheless, we 
will see that the added matrix structure (which encodes 
the chiral nature of the graphene electronic excitations) 
gives rise to qualitatively new effects. 



C. Diffusion propagator and renormalization to the 
e-phonon vertex 



1. Diffusion propagator 



It follows from Eqs. (11) and (21) that a key part of 



the disorder-induced modification of the electron-phonon 
heat flux is due to the modification of the effective 
electron-phonon interaction vertex. This modification is 
in turn directly related to the chiral diffusive dynamics 
of electrons in graphene, as described by Eq. (17). In 
this subsection, we discuss the diffusion propagator in 
more details, as well as the forms of the dressed electron- 
phonon vertices. Note that we restrict our discussion here 
(as we do throughout the paper) on electrons in the 
valley. While the sign of the chirality will be different 
for holes, or for the K_ valley, this sign has no impact 
on the quantity of interest, the e-phonon heat flux in the 
presence of disorder. 

We focus here on the most interesting diffusive regime, 
where ql ^ 1,ujt <C 1 (i.e., we are interested in length 
scales longer than / and time scales longer than r). In 
this limit we can work to lowest non-vanishing order in 
ql and ujT. The inverse diffusion propagator 'D~^(q,uj) — 



1 — 7'(q, uj) simplifies to 
il?-l(q,c.) = 
^ -iuj + Dq^ 



2 Vt 





iu - 





Dq^) 



2ivFqx 

2ivFqx D{ql-ql) 
2ivFqy ^Dq^qy 











2ivFqy 
2Dq^qy 



-Dq^ 









D{ql~ql) 







(23) 



in the diffusive limit, where D — vfI/2 is the usual dif- 
fusion constant in two dimensions. 

To gain intuition, it is useful to follow Ref. 1331 and con- 
sider the real-space representation of the matrix diffusion 
propagator, which describes the coarse-grained evolution 
of charge and pseudospin densities {N{r,t) and Sj{r,t) 
respectively, j = x,y,z). The first term in Eq. (23) 



would simply lead to uncoupled equations for each of 
these quantities: N would be described by a standard 
diffusion equation, while and Sy would have an addi- 
tional decay term (rate l/r), corresponding to the fact 
that pseudospin is not a conserved quantity. Sz has no 
dynamics in the limit we consider, as it precesses with 
frequency Ep and averages away on the timescale l/r. 
The second term in Eq. (23) complicates the above 



picture, as it now links the dynamics of charge and pseu- 
dospin densities. We thus obtain a set of coupled dif- 
fusion equations, describing the the dynamics of these 
quantities (note that we have taken into account the fact 
that the Pauli matrices are twice of the pseudo-spin ma- 
trices): 



= DV^N -vf{ 
3D 92 



9a;' 



dN 
~dt 

dt 2 
VF dN 
2" 9^' 
dSy ^ D(^^ 
dt 2 dx^ 
VF dN 
"T 9^- 



dSx 
dx 

D(P_ 
2 dy^ 



3p_cf_ 



dy 

Sx + D 



Sx + D 



92 

5* 

dxdy ^ 



Sx 



dxdy 



(24a) 

Sx 
T 

(24b) 

Sy 
T 

(24c) 



These equations are analogous (but not identical) to 
the diffusive dynamics for charge and spin in a 2D he- 
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lical metaP^ or the diffusion equations of Cooperons in 
graphend^. The interpretation here is similar to Ref . [Ml 
the coupling between charge and pseudospin dynamics 
in the diffusive limit is a result of the effective helicity of 
the electronic eigenstates. By helicity, we mean that at 
e > eigenstate of Eq. ([3|, pseudospin will be aligned 
with momentum. Thus, a positive gradient in say 
in the x-direction implies a corresponding positive gradi- 
ent in density of electronic a;-momentum. This will then 
naturally cause the charge density N to decrease in time: 
this is the third term in Eq. ( 24a ) . Alternatively, writing 



Eq. ( 24a ) in the form of a continuity equation, 



dN 
~dt 



= -V-J, 



(25) 



one sees that the charge current density J has the form 



J = -DVN + vp{S^x + Syij). 



(26) 



The first term is the usual diffusive current, while the 
second term corresponds to a "drift" current driven by 
the pseudospin density following from the Hamiltonian 
in Eq. (|3| (i.e. the current operator is the pseudospin 
operatorj! 

Turning to the dynamics of pseudospin densities. 



Eqs. ( 24b ) and ( 24c I again reflect the fact that pseudospin 



is not conserved, and effectively decays on a timescale r 
due to elastic impurity scattering. In addition, we see 
that the diffusion of these densities is anisotropic: this 
is also a simple consequence of helicity, as a net pseu- 
dospin density in a specific direction also implies a net 
momentum density in this direction which reinforces the 
diffusion in this direction. 



Finally, inverting Eq. (231 (using as always the Pauli 
matrix representation defined in Eq. ( |19[ )), one finds the 
diffusion propagator in Eq. ( 20 ) in the diffusive limit to 
be 



2?(q,c^) = 
1 

{-iu+2Dq^)T 

iqxl 3 
{iuj-2Dq^}T 2 

iqyl 
{iu}-2Dq^)T 





iq^l 
{iul-2Dq'^)T 
iujr+(ql-ql)l'' 
^ 2{iu)-2Dq'^)T 

q^qyi^ 

{iul-2Dq'^)T 





iqyl 



(iuj-2Dq^)T 

(iuj-2Dq'^)T 
ii^T-{ql-ql)l 
^ 2{iuj-2Dq^)r 

1/ 








2. Disorder vertex correction of deformation-potential 
e-phonon vertex 

Having discussed the basic form of the matrix diffu- 
sion propagator, we now turn to the renormalization of 
the e-phonon interaction vertex due to diffuson, as given 
by Eq. (191 and Eq. ( 27 ) . Consider first the DP contribu- 



tion to the vertex. The bare DP vertex is just a diagonal 
matrix in the sublatticc basis (c.f. Eqs. ([6])) proportional 
to the coupling constant gi. The impurity-dressed ver- 
sion takes the form: 



M, 



DP.diff 



(q,w) 



iqiq9i 



(~-iu+2Dq^) 
(-jw+2£)g2) {-iLu+2Dq^) 



iqle^^'*'i 
' {~iu:+2Dq^) 
1 



(28) 

The diagonal parts of the vertex acquire a diffusion pole 
protected by charge conservation, analogous to the case 
of a normal metal. We will be interested in Eq. ([28]) with 
w, q corresponding to a thermal phonon, huj = hsq ~ 
k^T. As s vp in graphene, we thus have that over a 
wide range of temperature 



VF 



-Tdis < T < Tdis, 



(29) 



one has uj ^ Dq^ , and thus the diagonal parts of the ver- 
tex in Eq. (28) will be enhanced by a factor ^ l/{q^P) 
compared to the clean case. This corresponds to the 
well known diffusive enhancement of the electron-phonon 
interaction: the diffusive motion effectively enhances 
the interaction time between an electron and a long- 
wavelength phonon. 

More surprisingly, Eq. ( 28 I implies that the diffusive 



renormalization of the deformation potential induces a 
vector potential which is along the direction of wave vec- 
tor q. This is a direct consequence of the chirality of the 
graphene electrons, which links the dynamics of charge 
and pseudospin, and thus allows a scalar potential to gen- 
erate a vector potential (i.e. a potential which couples to 
pseudospin). In the next section, we will show that this 
induced vector potential only gives a small contribution 
to the heat flux compared to the renormalized deforma- 
tion potential coupling. 



(27) 



Note that the effective diffusion constant (i.e. the coef- 
ficient of q^ in the diffusion poles appearing above) is 
twice the value of the standardly-defined D appearing in 
Eqs. (p4|): Doff = 2D = vpl- This effective doubling of 



the diffusion constant is a direct consequence of the chi- 
ral nature of electrons in graphene, and i s cons istent with 
the results of previous transport studie^^^'^. Also note 
that as expected from our discussion following Eq. ( 23 1 , 
only the charge-charge component of V (i.e. the (1, 1) 
matrix element) diverges in the limit of small lj and q. 
The lack of any corresponding large enhancement of the 
spin components of V is directly tied to the fact that 
pseudospin is not a conserved quantity. 



3. Disorder vertex correction of vector-potential e-phonon 
vertex 

The bare VP e-phonon vertex is given in Eqs. (|6|; it is 
purely off-diagonal in pseudospin space, and implies that 
a phonon of wavector q = gq = g(cos 0q, sin 0q) and po- 
larization 77 generates an effective vector potential Aq = 
52(sin20q,cos2(/)q) {jq = I) or A^ = g2(cos 20q, sin 20q) 
[r] = t). It is useful to decompose these vectors into their 
longitudinal and transverse parts: 



A '' - 4 '' n - 



(30) 
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FIG. 2: Color online. Thermal conductance per unit area 
G/A associated with the deformation potential coupling ver- 
sus temperature T, including the effects of disorder, but with- 
out electronic screening. We have taken a bare coupling con- 
stant g\ = 20 eV, carrier density n = 10^'^/cm^ and mean 
free path I = 20nm. The black solid line is the full result of 
our theory. The green-dashed line shows the asymptotic 
dependence in the low-temperature T <^ Tjis limit, whereas 
the red-dotted line shows the asymptotic behaviour in the 
high-temperature (clean) limit. 
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FIG. 3: Color online. Thermal conductance per unit area 
G/A associated with the deformation potential coupling ver- 
sus temperature T, including the effects of disorder, show- 
ing the effects of varying the mean-free path I as indicated; 
screening is neglected. Remaining parameters are the same as 
Fig. [2] Both the enhancement of the low-temperature ther- 
mal conductance and shift of the cross-over temperature with 
increasing disorder are clearly evident. 



where qj^ is a unit vector perpendicular to q. One finds 

^||,o = .92 sin3(/)q, A'-j_ q = 52 cos30q, (31a) 
^11 = .92 cos (/)q, A*j_ Q = -.92 sin 0q. (31b) 



Given the linearity of the vertex correction described 
by Eq. (17 1, we can separately analyze how disorder 



changes the interaction with the transverse and longi- 
tudinal phonon-induced vector potentials. Each of these 
will yield an e-phonon vertex which is a 2 x 2 matrix in 
pseudospin space. The bare vertices are: 



M^Q = (cOS0q(Ta; + sin^qfT;;) , (32a) 

-^1,0 ^ g (-sin(j!)q(7j; + cos(/)q(Tj,) . (32b) 



The corresponding renormalized vertices are given by 
Eqs. (Il9]) and ((27l), yielding: 



m; 



.diff 



2{iu: - Dq^) ^ 
{lL0-2Dq^) II 



M'' - 2M'' 



zZJ^^W"^"'"^""' ^^^^^ 
(33b) 



We see that the diffusive renormalization simply doubles 
the e-ph vertex associated with the transverse part of 
the phonon-induced vector potential; there is no diffu- 
sion pole here, as there is no charge associated with a 
transverse vector potential. In contrast, the longitudinal 
part acquires a diffusion pole. The vector potential part 
of the renormalized vertex is simply the bare vertex mul- 
tiplied by a factor of ^^^'^"^^2]^ ■ This factor tends to 
2 in the dc limit q — >■ 0, a; — >■ and 1 in the regime we 
are most interested in in this work, i.e. ql > s/vp. The 
renormalization factor of 2 in the dc limit is consistent 
with the re norma lization of current vertex in graphene in 
the dc liniii|2212Sl. Finally, we see that the vertex associ- 
ated with the longitudinal vector potential also acquires 
a scalar potential (second term in Eq. ( 33a I : this is an in- 
duced deformation potential coupling, again arising from 
the helicity of electrons in graphene. 

We stress that in contrast to the renormalized DP 



(c.f. Eq. (28)), Eqs. (33 1 explicitly show that there is 
no large enhancement of the VP vertex in the w = 
q — limit of interest. As discussed, the lack of a dif- 
fusive enhancement is a direct consequence of the non- 
conservation of pseudospin and the consequent lack of a 
protected diffusion pole. The net result is that the diffu- 
sive vertex correction discussed here does not significantly 
enhance the heat flux associated with the VP coupling 
at low temperatures. 

Finally, for completeness, we give the full form of the 
dressed VP e-phonon vertex. For the interaction with 



LA phonons, combining Eqs. (30 1 to (331 yields: 
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iql 

-iujT-\-q^l'^ 



sin 3(j)q 



y2 lUJT—q^l^ J —iujr-\-q^l^ 



r/2_ 



iql 



sin 30q 



(34) 



The full interaction vertex for TA phonons can be ob- 
tained in a similar fashion. One finds that vector po- 
tentials arising from TA and LA phonons make identical 
contributions to the heat flux. 



III. RESULTS AND DISCUSSIONS 



A. Heat flux without screening 



F„(T) {a = DP, VP), defined via: 



Having now determined both the 
electron-phonon vertices (c.f. Eqs. (28), 



renormalized 
(|34|) as well 



as the disorder-averaged electronic Green functions 
(c.f. Eq. (15l), we have all the necessary ingredients 



to evaluate Eq. (21) for the electronic collision integral. 
From this, Eq. (Ill directly yields the desired electron- 
phonon heat flux. As with the disorder-free case, we 
again find that the DP and VP couplings contribute in- 
dependently; we can thus meaningfully discuss the flux 
associated with each coupling. It is useful to express each 
of these heat fluxes in terms of an energy control function 



-Pa (^c 1 



ph) 



de£/a,o(e) = F„(re)-i^a(rph). (35) 



Here, Ia,o{£) is the momentum-average of the collision 
integral corresponding to the coupling mechanism a. 

We discuss each mechanism in turn, focusing as always 
on the regime ^Tdis <T <^ Tbg: as discussed, the lower 
limit here allows us to ignore the frequency dependence 
of the renormalized electron-phonon vertices. We also 
first discuss our results in the absence of any electronic 
screening. 



1. Deformation potential heat flux 

The energy control function determining the DP heat 
flux through the whole temperature regime ^Idis < 



T < Tbg is obtained from Eq. (jC^j) and Eq. (jCSj) to 
be 



vp 2ttpm J a 



dqq^ 



ql 



^/l + qH^ 



1 

ql 



1 - 



(36) 



Here, ujq = sq, and we have included an overall factor of 4 
reflecting the identical contribution from both valleys and 
both physical spin projections. The integrand reflects the 
contribution from phonons having momentum q to the 
heat flux. The first term corresponds to the contribution 
stemming only from the diagonal parts of the renormal- 
ized DP electron-phonon vertex, whereas the second term 
corresponds to off-diagonal terms (i.e. the effective vec- 
tor potential generated by the chiral diffusion) . The clean 
limit can easily be obtained by taking Z — >■ oo, yielding: 



DP, clean 



= lim i^Dp(T) = 



15^^ h^PMVpS 



3 „3 



(37) 

where we explicitly include factors of h in the expression. 
This result is the same as the heat flux for deformation 
potential in clean graphene obtained in previous theoret- 
ical work.-Siaiio] 



For finite I, both terms in Eq. (36) contribute. For 



long- wavelength phonons satisyfing ql 1, the first term 
dominates the second (effective vector potential) term by 
a large factor 1/(17/)^. In contrast, for phonons having 
ql ^ 1, both terms make comparable contributions. In 
the diffusive limit T < Tdis, the energy control function 
reduces to 



F m-2C(3)2 
-f'DP(J ) - — ^9i- 



(38) 



h'^PMs'^V%l 

The power law becomes instead of T"'; this is the result 
of the diffusive enhancement of the diagonal parts of the 
renormalized DP vertex, in complete analogy to what 
happens in a conventional disordered metal^^^^*^. The 
corresponding thermal conductance G = {d/dA)P{T + 
A, T) associated with the DP coupling is shown in Figs. [2] 
andlH 

As noted in the introduction, the expression in Eq. ( |38[ ) 
is exactly twice of that of the e-phonon heat flux as- 
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in a second-order process involving both e-phonon and 
e-impurity scattering events. 



sociated with impurity-assisted "supercolhsions" , as de- 
scribed by by Song et al. in Ref. [1] for temperatures 
T ^ Tbg • UnUke our study, Ref. [T] considers e-phonon 
scattering dressed by only a single impurity scattering 
event, and simply sums the effect of each impurity. In our 
approach, the g-dependence of the integrand in Eq. (36) 
in the dirtly limit I — can ultimately be traced to the 
diffusive enhancement of the DP vertex (c.f. Eq. (28 1). 
In contrast, the corresponding g-dependence in Ref. 1 

can be traced to the energy of a virtual electronic state Tbg is obtained from Eq. (|C4| and Eq. (|C5| to be 



2. Vector potential coupling heat flux 

The energy control function Fyp{T) for vector poten- 
tial coupling in the temperature regime j^Tdis < T <C 



FvpiT) 



452 



vf 27rpM 



dqq^ 



ql 



2ql 



(1 



v/i + g2p- 



iVK,T). 



(39) 



We have again included an overall factor of 4 in front to 
take into account the valley and spin degeneracy. The 
first term in the square bracket is the contribution from 
the vector potential after renormalization and the second 
term is the contribution from the induced deformation 
potential coupling from vector potential. At ql <^ 1, 
the second term gives a contribution about 1/6 of the 
first term, while at ql ^ 1, the second term is smaller 
by a factor of l/2ql. In the clean limit T ^ Tdis, the 
energy control function for vector potential has the same 
functional form as that for the DP coupling, i.e., 



Fvp{T) 



{kBTT 



(40) 



In the opposite diffusive limit T <C T^is, the energy con- 
trol function for the VP coupling reduces to 



^yp(J ) - „2 92 



EfI /, x5 

i6 [kBT) 



(41) 



Unlike the DP coupling, we see disorder increases the 
power of temperature of the low-temperature heat flux, 
indicating a suppression of heat flux from VP coupling. 
As discussed, the lack of a diffusive enhancement of the 
VP coupling is directly tied to the non-conservation of 
pseudospin; the main remaining effect of disorder is a 
simple broadening of the electronic eigenstates. The 
result is that for T <C T^is, disorder suppresses VP- 
mediated heat transport. Note that as the bare VP cou- 
pling constant 32 has been estimated to be more than 
an order-of-magnitude smaller than the corresponding 
DP coupling constant gi^^, it follows that in absence of 
screening, the heat flux associated with the VP coupling 
is expected to be negligible in comparison to that asso- 
ciated with the DP coupling. The thermal conductance 
associated with the VP coupling is shown in Figs. |4]and 

m 



B. Heat flux with electronic screening 

We now consider how the above results are altered if 
one includes the screening of the e-phonon interaction. 
As discussed extensively by von Oppen et aP^, the defor- 
mation potential coupling will be subject to screening at 
long wavelengths in the usual manner, whereas the vec- 
tor potential will not be screened, as it does not induce 
any net charge density. In this work, we are interested 
in temperatures such that T > (s/wF)7dis (c.f. Eq. ([29|)), 
implying that dynamic screening effects (which are also 
sensitive to disordeil^ will be unimportant. Further, in 
the temperature regime of interest (T < Tbg), a sim- 
ple Thomas-Fermi approach to screening is expected to 
suffic^2ll_ The result is that the bare DP coupling con- 
stant gi in Eqs. ^ now becomes dependent on the mag- 
nitude of the phonon wavevector q: 



5i,sc(g) = 51- 



— , (42) 
q + qxF 

where the Thomas-Fermi wavector gxF is given bjE^ 



(7tf 



KhVF 



-kf 



(43) 



and K is an effective dielectric constant. Using the value 
of K appropriate to graphene on a Si02 substrate, one 
has qtf — 3.2fci7^^ Note that as we focus on the regime 
I^fI ^ 1, the effects of screening will generally set in 
at a much higher temperature Tgc = s^tf/^b than the 
temperature Tdis below which disorder-effects becomes 
important. We note that a recent experiment measuring 
the e-phonon contribution to the electrical resistivity of 
a suspended graphene flake suggests that screening does 



not seem to be playing a role even when T < Ts 



1261 



as 



the results are compatible with the predictions for an 
unscreened deformation potential interaction (see also 
Ref.l 
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FIG. 4: Color online. Thermal conductance per unit area 
G/A associated with the vector potential coupling versus tem- 
perature T, including the effects of disorder, but without elec- 
tronic screening. We have taken a bare coupling constant 
g2 = 1.5 eV, carrier density n = 10^^/cm^ and mean free 
path I = 20nm. The black solid line is the full result of our 
theory. The green-dashed line shows the asymptotic de- 
pendence in the low-temperature T <^ Tdia limit, whereas 
the red-dotted line shows the asymptotic behaviour in the 
high-temperature (clean) limit. 



One can now easily include the effects of screening 
into our theory by making the substitution gi — >■ gi,sc{q) 
in Eq. (36) for the energy control function F£,p{T) de- 
termining the deformation-potential mediated heat flux. 
One finds: 



DP, I 



:(r) = 



_ I 63p, 



24 gfEf 



T6, 

C(5)r5, 



if Tdis « T « Ts, 
if T < Tdis,Tsc- 



(44) 

The suppression of the DP coupling by screening at 
low temperatures T <C T^is implies that its associated 
heat flux can now become comparable or even smaller 
in magnitude to that associated with the VP coupling, 
c.f Eq. (41). In the low temperature limit, both mecha- 
nisms yield energy control functions F{T) cx T^, with: 



i^DP,sc(T) gl 1 



Fyp{T) gl qlpP 



(45) 



We see that the largeness of gxF^ can compensate for 
the relative smallness of §2 with respect to 171, lead- 
ing both mechanisms to make comparable contributions. 
This behaviour is demonstrated in Fig. [6j where the 
thermal conductance versus temperature for both mech- 
anisms is presented, for both the strongly and weakly- 
screened cases. The fact that both mechanisms are 
comparable is markedly different from what happens in 
the screened, disorder-free case, which is realized when 
Tdis < T < Tsc. In this case, the VP heat flux wifl dom- 
inate the DP heat flux at low temperatures, as it scales 
like as opposed to (see Fig. [7|. This behaviour 
in the clean limit is similar to expectations for e-phonon 



FIG. 5: Color online. Thermal conductance per unit area 
G/A associated with the vector potential coupling versus tem- 
perature T, including the effects of disorder, showing the ef- 
fects of varying the mean-free path I as indicated; screening 
is neglected. Remaining parameters are the same as Fig. [4] 
Both the suppression of the low-temperature thermal conduc- 
tance and shift of the cross-over temperature with increasing 
disorder are clearly evident. 



contribution to the electrical resistivity, where it has also 
been argued that the VP coupling can dominate at low 
temperature^. 

Finally, we note that with our theory, the only way 
to obtain a power law in the heat flux at low tem- 
peratures (i.e. ^ = 3 in Eq. ([T])) is via an unscreened 
deformation potential. Thus, measurements of the low 
temperature heat flux could also serve as a diagnostic 
tool for assessing the importance of screening the defor- 
mation potential. 



IV. CONCLUSIONS 

We have presented a comprehensive theory showing 
how electronic disorder modifies the electron-phonon in- 
teraction in graphene (both the vector potential and de- 
formation potential couplings), and how this in turn has 
observable consequences for the heat flux between the 
electrons and lattice (acoustic) phonons. We focused 
on the relatively simple situation where the graphene is 
doped away from the Dirac point, and where the impu- 
rity potential can be considered smooth on atomic scales 
(implying that the disorder potential preserves the pseu- 
dospin and valley symmetries of the graphene Hamilto- 
nian). We found that the unusual diffusion dynamics 
of electrons in graphene that results from their chirality 
also has implications for how disorder modifies electron- 
phonon physics. We also found that this modification is 
quite different for the deformation potential coupling ver- 
sus the effective vector potential coupling. In the absence 
of screening, the contribution to the heat flux in Eq. (fTl) 
from both couplings has temperature dependence of 
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FIG. 6: Color online. Thermal conductance per unit area 
G/A versus temperature T, showing the effects of electronic 
screening, gi — 20 eV, g2 = 1.5 eV, a carrier density 
n = Ifli^^/cm^ and a mean free path I = 20 nm. The red 
short-dashed line correspond to an unscreened deformation- 
potential coupling, while the blue dashed-dotted line corre- 
sponds to a screened deformation-potential coupling, with a 
Thomas-Fermi wavevector qtf = 3.2kF as appropriate for 
graphene on SiOf^. The dashed purple curve is the contri- 
bution from the vector potential coupling. Despite its much 
smaller bare coupling constant, we see that at low temper- 
atures, both the deformation potential and vector potential 
couplings make almost equal contributions when both screen- 
ing and disorder effects are included. 



in the clean limit T 3> s/l, consistent with previous work. 
In the disorder limit T < s/l, however, disorder affects 
the two types of couplings differently and the power law 
of heat flux for the two couplings becomes different. We 
found that the total effect of disorder enhances the heat 
flux from DP coupling, however, it suppresses the heat 
flux from VP coupling. The power law in Eq. ([I]) for 
DP coupling becomes in the disorder limit while for 
VP coupling, the power law becomes in the same limit 
without screening. Thus, without screening, the DP cou- 
pling is expected to dominate heat transport in both the 
clean limit and disorder limit, given that DP coupling is 
believed to be much stronger than vector potential cou- 
pling. 

We also considered the effects of screening on the 
above physics, adopting the Thomas-Fermi approxima- 
tion. Only the DP is affected by screening; it suppresses 
it, and thus at low temperatures, its contribution to the 
heat flux scales like in the clean limit and in the 
disordered limit. Without disorder, one would thus ex- 
pect the VP to dominate at low temperatures due to the 
screening of DP, similar to expectations for the electron- 
phonon contribution to the resistivity in clean graphene 
sheets^'^. However, when disorder effects are included 
along with screeening, we flnd that both the DP and VP 
coupling mechanisms can make comparable contributions 
at low temperature. 

In the future, it would be interesting (though challeng- 
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FIG. 7: Color online. Thermal conductance per unit area 
G/A versus temperature T, showing the effects of electronic 
screening. Parameters are identical to Fig. [6] except we have 
now taken the mean free path to be Z = 10 ^m, meaning that 
we are effectively in the clean limit. Unlike the disordered 
case shown in Fig. |6] we now see that at low temperatures, 
the contribution of the vector potential coupling dominates 
that from the deformation potential coupling. 



ing) to extend these results to situations of lower doping 
(where the Fermi energy is closer to the Diract point), 
and to include a richer structure of disorder that can 
break the symmetries of the clean graphene Hamiltonian 
(as has been done in, e.g. , studies of weak-localization-'-^') . 
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Appendix A: Keldysh formalism of the kinetic 
equation 



The Green's function in the Keldysh space is a matrix 



as 



G 




(Al) 



where and G^ are the retarded and advanced Green's 
function respectively. The Keldysh component G^ in 
general could be parameterized as G^ = G^oF — FoG^ , 
where F is the Hermition distribution function matrix. 
The Green function obeys the following Dyson's equa- 
tion: 



(G^)-i - 
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where the bare retarded and advanced Green function in 
graphene are (Gq^^)~^{x, x') = S{x — x'){idt> + id ■ Br') 
in the sublattice basis. Yj^/^ and are the retarded, 
advanced and Keldysh component of the electron self en- 
ergy respectively. The circle o indicates integration over 
intermediate coordinates. 

The equation for the Keldysh component reads 

Fo [G^]-^ - {G^y^ oF + Fo±^~±^oF. (A3) 

The distribution function matrix F{x,x') in graphene 
could be decomposed to four components, 



(A4) 



i=0 



where a'^ is the 2x2 unit matrix and (7*,i — 1,2,3 are 
the Pauli matrices. 



The Wigner transformation of Eq. ( A3 ) gives the ki- 



netic equation of the distribution function as: 

3 3 

-'^ia'drfi{t,p;e,p) -2ivF ^ e^jkd''kjfi{t, p;e,p) 

i—G i,j,k—l 

= p) + Fit, p- e, p)S^(e, p) - p)F(t, p- e, p) 

= /(e,p;t), (A5) 

where the function fo{t, p; £,p) = 1 — 2n(i,p;e,p) and 
n(t, p;e,p) is the charge density distribution function, 
fi(t, p;e,p), i — 1,2,3 represent the pseudo-spin den- 
sity distribution and e^j-fc is the three-dimensional anti- 
symmetric tensor. The parameters t and p are center of 
mass time and coordinates respectively. Since the trans- 
lation symmetry in space is restored after averaging over 
impurities, we drop the t dependence in the text. 

The right hand side of this equation is the collision 
integral I{e, p; t) (times i) in presence of interactions. 
T,^- ,1]^ ,T,^ are respectively the retarded, advanced and 
Keldysh components of the electron self energy due to 
e-phonon and impurity scatterings. From the left hand 
side of the above equation, one finds that the pseudo-spin 
density distribution functions /i=i,2,3 are small in a fac- 
tor of l/EpT compared to fo- To leading order, one can 
replace F{t,p;e,p) on the right hand side by the scalar 
charge component fo{e,p;t). 

The collision integral is a 2 x 2 matrix in the sublattice 
basis, which can be decomposed to components of (7° and 
CT* [i = 1, 2, 3). The three components of ct* (i = 1, 2, 3) 
give a measure of the pseudo-spin density induced by 
interaction, while the a'^ component determines the col- 
lision integral for the charge density distribution function 
and is the one of interest in this paper. From Eq.(A5l, 



one gets the collision integral for charge distribution func- 
tion 



1 9/o(e,p;t) 

2 dt 



-|Tr 



(£, p) + (1 - 2n(£, p; t))(S^ - E«)(e, p) 



(A6) 



Appendix B: Full form of the diffusion propagator 
and renormalized e-phonon vertex 

The full form of the diffusion propagator X'(w, q) cross- 
ing the whole temperature regime ^Jdis < T ^ Tbg is 
quite complicated, yet in the temperature regime T > 
:^Tdis, the frequency dependence of the diffusion prop- 
agator can be dropped and the diffusion propagator in 
Eq. ( 20 ) is simplified to 



as shown in Eq. ( 12 I 



>0,q) = 

-^COS0q 

-^sin^iq 



COS 

ql 
1 



q 

cos 20q 



2^1+92(2 
sin 20q 





-^sm(^q 

sin 20q 
2^1+92;2 

, 1 + COS 2</)q 










1/ 



(Bl) 

in the Pauli matrix basis. 

The four vector representation of the bare deformation 
potential vertex in the Pauli matrix basis is mop, o = 
jg^' (1, 0, 0, 0)(7i. The renormalized vertex according to 



Eq. ( 19 ) in such basis becomes 

(l 

m-DF, diff(w 0, q) = iqS} 



r COS (br, 



V 



sine 




51, (B2) 



where the superscript T means transpose of the column 
vector to row vector and the same for myp, diff below. 
Written in the sublattice basis, the renormalized defor- 
mation potential vertex becomes 



MDP^diff(w 



►0,q) 

q^l^ 

ie''t'<i 
ql 



1 



ql 

q-^l-^ 



91- (B3) 



In the diffusive limit lot <^ l,ql <^ 1, it reduces to 
Eq. ( 28 ) (dropping the frequency dependence there) ; 
while in the clean limit ^ 1, it reduces to the bare 
vertex. 

The four vector representation of the bare vector po- 
tential for LA phonon in Eq. Q is 

mvp,o = iqS.1 {0,sm2(j)q,cos2(l)q,0) g2 (B4) 
in the Pauli matrix basis. The renormalized vertex ac- 



cording to Eq. ( 19 1 then becomes 
TOvp,diff(t^ 0, q) 
/ 



sin 2(/)q(l 

cos 2(/)q(l 



i sill 30q 

ql 

1 



r 



1 



)- 



V 



sin 4(/)q 
COS 4(;6q 





92 



(B5) 
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in the same basis. Written in the sublattice basis, it reads 



Mvp,difr(t^ 0,q) 



i sin 3^q 



i sin 30q 



(B6) 



In the diffusive limit, it reduces to Eq. (34) (again drop- renormahzed deformation potential to Eq. (21) and inte- 



ping the frequency dependence there); while in the clean grating over the phonon frequency, one gets the collision 
limit ^ 1, it reduces to the bare vertex. integral for deformation potential coupling as 



Appendix C: Details of the matrix collision integral 

a. Deformation potential 

The renormahzed deformation potential in the whole 
regime of ql is presented in Appendix B. Plugging in the 



^DP,o(e,p) 



dn{e,p;t) . f dq 



dt 



(2vr) 



(e + t^g+27)'-4lP + ql' 



(CI) 



where LOq = sq. 

The collision integral for deformation potential after 
average over the electron momentum becomes 



-fDP,o(e) = / 7;^A(£,p)/dp,o(£,p) 



■KV J (27r) 



' (1- . ^ ) 



q^P ' ^/l + qH^ q^P' ^/l + Wl 



i?(e, ujq) - {ujq -ujq, q -q)}. 



(C2) 



The heat flux is then 



where the energy control function Fdp{T) is presented 
in Sec. UHl 



PDpiTc,Tpi,) = 1/ / dee/DP,o(£) = FnpiTp^) - Fdp{T, 



(C3) 
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6. Vector potential the component for the charge distribution function, one 

gets the kinetic equation of the charge distribution func- 
Plugging in the renormahzed vector potential in Ap- tion due to vector potential coupling after average over 
pendix B to the collision integral Eq. ^ and separating the angel of electron momentum as 



ot J (27r)^ 2^1 + ^2/2 + 2q^l'^ ^i^qip 

-(w^ -tj„q^ -q)}. (C4) 



The heat flux due to vector potential coupling is where the energy control function i^vp {T) is presented in 

Secinil 

PvpiT,,Tpi,) j d£e/vP,o(e) - i^vp(re) - Fvp(Tph), 

(C5) 
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